arXiv:1506.08719vl [astro-ph.EP] 29Jun2015 


Draft version June 30, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


FOLLOW-UP OBSERVATIONS OF PTFO 8-8695: A 3 MYR OLD T-TAURI STAR 
HOSTING A JUPITER-MASS PLANETARY CANDIDATE 

David R. CiardP, Julian C. van Eyken^’^’^, Jason W. Barnes'^, Charles A. Beichman\ Sean J. Carey®, 
Christopher J. Crockett®, Jason Eastman®’’’, Christopher M. Johns-Krull®, Steve B. Howell®, Stephen 
R. Kane’®, Jacob N. Mclane’I Peter Plavchan’®, L. Prato”, John Stauffer®, Gerard T. van Belle”, Kaspar 

VON Braun” 

Draft version June 30, 2015 

ABSTRACT 

We present Spitzer 4.5/im light curve observations, Keck NIRSPEC radial velocity observations, 
and LCOGT optical light curve observations of PTFO 8-86 95, which may host a Ju pite r-sized planet 
in a ve ry short orbital period (0.45 days). Previous work bv Ivan Evken et 'HI (j2012h and iBarnes et al.l 
(|2013f) predicts that the stellar rotation axis and the planetary orbital plane should process with a 
period of 300 — 600 days. As a consequence, the observed transits should change shape and depth, 
disappear, and reappear with the precession. Our observations indicate the long-term presence of the 
transit events (> 3 years), and that the transits indeed do change depth, disappear and reappear. The 
Spitzer observations and the NIRSPEC radial velocity observations (with contemporaneous LCOGT 
optical light curve data) are consistent with the predicted transit times and depths for the M* = 

0.34 Mq precession model and demonstrate the disappearance of the transits. An LCOGT optical light 
curve shows that the transits do reappear approximately 1 year later. The observed transits occur at 
the times predicted by a straight-forward propagation of the transit ephemeris. The precession model 
correctly predicts the depth and time of the Spitzer transit and the lack of a transit at the time of 
the NIRSPEC radial velocity observations. However, the precession model predicts the return of the 
transits approximately I month later than observed by LCOGT. Overall, the data are suggestive that 
the planetary interpretation of the observed transit events may indeed be correct, but the precession 
model and data are currently insufficient to confirm firmly the planetary status of PTEO 8-8695b. 

Subject headings: (stars:) planetary systems - stars: individual (PTFO 8-8695, 

2MASS J05250755-f0134243, CVSO 30) - stars: pre-main sequence 


1. INTRODUCTION 

With the discove ries fueled by the Kepler Mission 
(|Borucki et al.l[20TOl . there are now more than 1800 con¬ 
firmed or validat ed planets (e.g., iBatalha et al.l I201II : 
iRowe et al.l 120141 1. Kepler has increased our knowl¬ 
edge of the diversity of planets and planetary systems, 
and the sheer number of planets discovered by Ke¬ 
pler, coupled with on-going discoveries from other tran¬ 
sit programs, radial velocity surveys, dir ect imaging 
effort s , and microlensing campaigns (e.g.. iGillon et al.l 
2014 Iwittenmver et al.l [20141 : [K rails fc Ireland 120121 : 
Koshimoto et al .I I2014I1 . tiave spawned a realization that 
nature yiel ds planetary systems in a variety of configura¬ 
tions (e.g..lWinn fc Fabrvckvll20l4lFabrvckv et al.ll20l4 
iCiardi et al.l 120131 : iSteffen et al.l 1201311 . To understand 
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the processes that shape the planetary systems that we 
observe today, it is crucial to observe the planets during 
the period of formation and evolution. 

Evidence has been found for extra-solar rings around 
a 16 Myr-old pre-main sequence star in the Upper 
Centaurus-Lupus subgroup of Sco-Cen, possibly indica- 
tive of early planet f ormation around a T Tauri star 
(IMamaiek et al.l [201211 . However, discovering and ob¬ 
serving planets in the earliest stages of formation and 
evolution requires observing stars that are only < 5 — 
10 Myr old, as this is the timescale ov er which the 
planet forming disks a re depleted (e.g., iHaisch et al.l 
120011 : iHillenbrandl l2008l l . although recent work may in- 
dicate that the disk s may last as long as < 20 Myr 
(|Pfalzner et al.l [201411 . But stars at this age are notori¬ 
ously active and spotted, making them photometrically 
and spectroscopically variable and ma king planet discov¬ 
ery at these young ages difficult (e.g., iMiller et al.l 120081 : 
iCrockett et al.l[2^^ iCodv et al.ll2014ll . 

The Palomar Transient Eactory (PTF) Orion sur- 

vey, conducted in _ th e win ters of 2009 and 2010 

(Ivan Evken et al.l 120111 1201211 . was an attempt to ad¬ 
dress this deficit by searching for young transiting hot 
Jupiters in the 25 Ori g roup, a small association o f T- 
Tauri stars ident i fied b y iBriceho et al.l (120(1511200711 . In 
Ivan Evken et al.l (|2012ll . we reported the discovery of a 
promising young planetary candidate orbiting the known 
weak-lined T-Tauri star PTF O 8-8695 with an ag e pre¬ 
viously estimated of « 3 Myr (|Bricefio et al.l [200511 . Su- 
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perimposed on top of notable quasi-periodic variability, 
the star showed regular transit events with a period of 
just 0.45 days. But unlike TW Hy a, which is a classical 
T Ta uri star with active accretion (jRucinski fc Krautterl 
llOS,*!!) and has periodic ~ 2% dips likely cau sed by the 
disk occulting a hot spot (jSiwak et al.ll20l4i . PTFO 8- 
8695, being a WTTS, has no infrared excess and likely 
no disk to occult or veil the star. Finally, no evidence 
was found for a stellar companion. 

The observed optical radial velocity variations (Keck 
HIRES and HET HRS) were found to be out of phase 
with the transit events and were not of sufficient ampli¬ 
tude to be caused by a stellar companion, arguing that 
the obs erved transits were not the result of an eclipsing 
binary (|van Evken et al.ll2012D . This is in contrast to 
V471 Tau which shows similar depth eclipses, but also 
150 km/s amplitude radial velocity variations in phase 
with the eclipses ca used by a stellar mas s white dwarf 
orbiting a K-dwarf (|Kamihski et al.ll2007f) . The radial 
velocity variations observed in the optical for PTFO 8- 
8695 are likely caused by the rotation of the spotted 
stellar surface (iPrato et ahl 120081: iHuelamo et al.l U008l 
iHuerta et al.1 I2008D . The spot induced observed ra¬ 
dial velocity variations in the optical ('^ 2 km/s semi¬ 
amplitude) yield a candidate companion mass upper 
limit of Msini < 4.8 Mj^p- Coupled with the transit 
modeling and the constraints on the line of sight inclina¬ 
tion {i Ki 60°. Ivan Evken et al.ll201^ . the upper limit of 
the companion mass was found to be Mp < 5 Mjup - 
well within the planetary mass regime. 

The transit light curves did show shape changes be¬ 
tween the two y ears in which the P TF Orion observa¬ 
tions were made. iBarnes et al.1 (|2013li showed that these 
effects could result from planetary transits of a fast- 
rotating, oblate, and significantly gravitationally dark¬ 
ened star, which is consistent with the short rotation 
period and high projecte d rotational velocity of P TFO 8- 
8685 (i;sini Ri 80 km/s. Ivan Evken et al.l[201^ . Being 
non-s pherical, and da rker at the equator than at the 
poles (Ivon ZeiDell924[i , oblique transits of such a rapidly 
rotatin g star can sho w unusual and asymmetric transit 
shapes ( BarneslF2009[ ) . Furthermore, if the planet orbit is 
oblique to the stellar rotation axis, a torque is exerted on 
the orbit by the stellar oblateness, leading to precession 
of the orbital nodes. 

Orbital precession has been noted previously in other 
systems and its effect on observed planet ary transits - 
parti c ularly in circumbinary s ystems (^e.g.. iLeung fc 
120131 : iMartin fc TriaudI I2014D . Recently, such behav¬ 
ior has p reviously be e n not e d in t h e Kep l er 13b planet 
syste m (|Szab6 et al.l I2011L 120121 12014 IBarnes et al.1 
1201 If), and with th e circumbinary planet, Kepler 413b 
(jKostov et al.l [ 201 ^ . These are predicted to proce ss on 
timescales of decad es. But Ivan Evken et al.1 (j2012f) and 
IBarnes et al.l (l2013[l showed that, for PTFO 8-8695, the 
precession period was on the order of a 300 — 600 days - 
a timescale acc e ssible to observation. 

IBarnes et al.1 (|2013[ ) were able to simultaneously fit 
both years of the PTF Orion data set using a fully 
consistent model which included gravitational darkening, 
oblateness of the host star, and an analytical model of 
the precession of the system. The fits yielded a planet 
mass Mp ^ 3 Mj^p, and predicted rapid precession on 
approximately year-long timescales. The fits implied a 


very high obliquity for the orbit (« 70°), resulting in a 
prediction that the transits should change depth, shape, 
and disappear for periods of a few months to a year as 
the precession brings the planet’s orbit out of the line of 
sight to the star. The limited timespan of the data and 
the uncertainty in the stellar mass and radius left an un¬ 
resolved degeneracy in the solution. More observations 
were needed to test the veracity of the model and better 
determine the system parameters. In particular, observa¬ 
tions were needed to test that the transits changed shape, 
disappeared, and reappeared in a predictable manner. 

Towards this end, we obtained follow-up observations 
of PTFO 8-8695 which included Spitzer photometry to 
confirm the transit events in the infrared; Keck NIR- 
SPEC infrared radial velocity measurements to mea¬ 
sure the radial velocity sign ature of the planet and the 
Rossi t er-McLaughlin effe ct 


11924 iQueloz et al.l 1200(11 : lOaudi fc WinnI 12007 ) as the 


iMcLaughlinI 


planet transited the star , hinted at in the orig inal op¬ 
tical radial velocity data (Ivan Evken et al.ll2012 ): optical 
photometric monitoring with the Las Cumbres Observa¬ 
tory Global telescope Network (LCOGT) to support the 
radial velocity measurements and to better establish the 
transit parameters and orbit. 

The Spitzer data, the Keck NIRSPEG data, and early 
optical LGOGT photometry data were acquired before 
the realization that the transits could disappear com¬ 
pletely, but these data, in the end, provided evidence for 
the changing and disappearing transits. Further optical 
LGOGT observations were obtained with the realization 
that the system could precess and that the transits could 
disappear. The optical data show the return of the tran¬ 
sits at the time predicted by propagating the transit mid¬ 
point ephemeris, but approximately 1 month prior to the 
prediction of the precession mo del. Given the complex - 
ity of the system as des cribed in Ivan Evken et aH (|2012h 
and IBarnes et al.l (120131) . we regard the partial matching 
of the observations to the model as an indication that 
the planetary interpretation of the nature of the source 
of the transits is still viable, but we are unable to confirm 
the planet with these data and the current models. 

2. OBSERVATIONS AND ANALYSIS 
2.1. Spitzer 4-5pLm 

Spitzer Space Telescope Director’s discretionary time 
was granted to observe PTFO 8-8695 on 28 April 2012 
UT, soon after the discovery paper was published. The 
primary purpose of the observations was to confirm 
the presence of the tra nsit two years after t he PTF 
Orion data were acquired (|van Evken et al.ll2012f) . and to 
search for a secondary eclipse. The infrared observations 
of PTFO 8-8695 consist of an 11.8 hour stare at 4.5^m 
with IRAG (request key 45476608), consisting of 1368 
thirty second integrations with an effective cadence of ~ 
31.3 seconds. The data were collected using the standard 
techniques for high-precision photo metry with IRAG t o 
minimize instrumental systematics (jlngalls et al.ll20O) . 
The observations were placed on the part of an array 
pixel that has minimal response variations. The stan¬ 
dard Spitzer Science Genter data products were used in 
our analysis. 

Aperture photometry was performed on each basic cal¬ 
ibrated data (BCD) image using a three pixel radius 
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TABLE 1 

Spitzer Light Curve Data (sample) 


BJD 

Phase 

Flux 

Uncertainty 

Flux Model 

Flux—Model 

Transit Model 

2456045.70197 

-0.0014 

0.9912959 

0.0032780 

0.9981057 

-0.0068097 

-0.0058880 

2456045.70233 

-0.0006 

0.9932295 

0.0033324 

0.9980942 

-0.0048646 

-0.0058900 

2456045.70269 

0.0002 

0.9914524 

0.0037372 

0.9980825 

-0.0066302 

-0.0058890 

2456045.70305 

0.0010 

0.9927249 

0.0039651 

0.9980709 

-0.0053461 

-0.0058860 
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Fig. 1. — Top: April 2012 Spitzer 4.5/rm light curve pl otted as a function of t ime. The blue circles indicate the data points predicted to 
be within the transit window from the PTF ephemeris llvan Evken et ani2012lh the red squares are outside the transit window and were 
used to fit an 8*^-order polynomial (black line) to model the stellar variability. Bottom: The Spitzer 4.5/rm light curve with the polynomial 
fit subtracted to remove the stellar variability is plotted as function of orbital pha se. The solid black line represents the transit fit to the 
light curve and the data have been phased on this new ephemeris (see Section [3T}. 


aperture and subtracting a background using an esti¬ 
mate of the mode of a circular annulus with inner radius 
of 3 and outer radius of 7. Source positions were de¬ 
termined on each BCD by calculating the center-of-light 
using a 5 X 5 pixel box centered on the peak pixel as¬ 
sociated with the source using the Spitzer Science Cen¬ 
ter provided IDL procedure box_centroid.pro; aperture 
photometry was performed utilizing the IDL version of 
the DAOPHOT routine aper.pro. Photometric system- 
atics were removed by applying an intra-pixel gain map 
([Ingalls et al.ll201^ to each photometric point as a func¬ 
tion of centroid position. The Spitzer light curve is tab¬ 
ulated in Table [Tj and shown in Figure [11 

2.2. Keck-II NIRSPEC Radial Velocities 

Keck-II NIRSPEC observations were obtained on 21 
December 2012 - a full two years after the detection of 
the optical transits with PTF. The purpose of the NIR¬ 
SPEC observations was to obtain phase-resolved high 
precision radial velocity data, primarily during the tran- 


TABLE 2 

Keck NIRSPEC Radial Velocity Data 


HJD 

Phase^ 

Radial Velocity 
[km/s] 

Uncertainty 

[km/s] 

2456282.791 

-0.272 

-0.971 

1.034 

2456282.805 

-0.241 

-0.098 

1.082 

2456282.818 

-0.212 

0.265 

0.932 

2456282.832 

-0.181 

0.364 

1.000 

2456282.872 

-0.092 

0.963 

0.975 

2456282.886 

-0.060 

0.890 

0.902 

2456282.902 

-0.024 

-1.104 

0.945 

2456282.916 

0.007 

1.398 

0.881 

2456282.934 

0.047 

-0.642 

1.065 

2456282.948 

0.078 

-0.456 

2.327 

2456282.964 

0.114 

0.800 

0.999 

2456282.977 

0.143 

-0.114 

1.046 

2456282.995 

0.183 

-0.866 

1.055 

2456283.009 

0.214 

-0.743 

1.798 

2456283.024 

0.247 

-2.783 

1.828 

2456283.038 

0.279 

-1.003 

1.167 


Based upon Spitzer Ephemeris (Table [4]l. 
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Fig. 2.— Top: December 2012 Keck NIRSPEC radial velocity data plotted as a function of time. The horizontal dashed line represents 
the weighted mean of the relative velocities and the solid black line represents the best fit radial velocity curve with the transit midpoint 
(vertical dashed line) fixed by the Spitzer ephemeris (Section 13.111 : the horizontal dotted line represents the velocity offset of the fit 
(7 = 0.237 ±0.415 km/s). The vertical dotted lines represent the predicted beginning and ending times of ingress and egress, respectively. 
Bottom: The radial velocity curve is phased on the Spitzer ephemeris. The overplotted lines represent the fits as in the top panel; the gray 
points are the values repeated in phase for continuity and clarity. 


sit of PTFO 8-8685b. The Rossiter-McLaughlin effect 
was predicted to be w 2 — 3 km/ s given the fsinz of th e 
star and the depth of the transit (jvan Evken et al.ll2012h . 

Spectra were acquired through the 0.432" (3 pixel) 
wide slit; the echelle and cross disperser angles were ori¬ 
ented to obtain K-band spectra containing Na I, Mgl and 
CO. In this configuration, orders 32-38 (non-contiguous) 
were imaged through the NIRSPEC-7 blocking filter. 
The data were acquired in pairs of exposures of 600 sec¬ 
onds in a standard ABBA sequence, with the telescope 
being nodded ±6" from the slit center along the 24" slit, 
so that each frame pair would contain object and sky in 
both nod positions. 

To correct for telluric absorption, we also obtained 
spectra of HR 2315, an AOV star located close to 
PTFO 8-8695 in the sky, before and after the PTFO 8- 
8685 observations. After every two ABBA sets, spectra 
of the internal NIRSPEC continuum lamp were taken 
for flat fields at the K-band settings; lamp exposures of 
the argon, neon, krypton, and xenon arc lamps provided 
wavelength calibration for all the K-band orders. At the 
end of the PTFO 8-8695 science observations, GJ281, a 
radial velocity standard, was observed to set the absolute 
velocity scale of the observations. The spectral images 
were dark subtracted and flat fielded and the individ¬ 
ual nods were extracted, wavelength calibrated with the 
lamp spectra, and telluric divided. Four ABBA nod sets 
were obtained for a total of 16 spectra and radial veloc¬ 
ity estimates. The night was plagued with variable high 


cirrus and as a result, the signal-to-noise per resolution 
element varied but was approximately 50. 

The Na I lines at ^ 2.206^m are the strongest lines 
in the near-infrared spectrum of PTFO 8-8695 and the 
entire order containing these lines was used to determine 
the relative radial velocities. The telluric corrected indi¬ 
vidual spectra were coadded to produce a master spec¬ 
trum with S/N ^ 200. The master spectrum was cross 
correlated with each of the 16 individual spectra to deter¬ 
mine the relative radial velocities between the spectra. A 
Monte Carlo effort was employed to estimate the uncer¬ 
tainties. The Monte Carlo simulations assume that the 
wavelength fits are perfect, which, of course, may not 
be correct. To try to characterize the wavelength fitting 
uncertainty, the spectra with strong telluric lines were 
utilized. The rms scatter in the radial velocity measure¬ 
ments of the telluric lines in each of the 16 spectra was 
found to be 0.68 km/s. This uncertainty contribution 
was added in quadrature to the Monte Carlo uncertainty 
to arrive at a final uncertainty. The measured radial 
velocities and associated uncertainties are presented in 
Table O and the radial velocity curve is shown in Fig¬ 
ure [2l 

2.3. LCOGT Photometry 

We obtained optical photometry for PTFO 8-8695 us¬ 
ing the LCOGT Im telescope network during Decem¬ 
ber 2012 and November 2013. We utilized the full net¬ 
work capabilities of LCOGT including the sites at Cerro 
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BJD - 2456284.0 [days] 



Fig. 3.— Top: December 2012 LCOGT r' light curve, obtained two days after the Keck NIRSPEC observations, is plotted as a function 
of time. The blue circles indicate the data points predicted to be within the transit window from the Spitzer ephemeris (Table : the 
red squares are outside the transit window and were used to fit a cubic spline (black line) to model the stellar variability. Bottom: The 
LCOGT r' light curve, phased upon the Spitzer ephemeris and with the spline fit subtracted, shows no sign of a transit. The solid black 
line represents a constant value of 0.0. 


TABLE 3 

LCOGT Light Curve Data (sample) 


BJD 

Phase 

Flux 

Uncertainty 

Flux Model 

Flux—Model 

Transit Model 

2012 December 

2456284.54447 

-0.4014 

1.0143089 

0.0052217 

1.0215716 

-0.0072627 


2456284.54568 

-0.3987 

1.0232391 

0.0051967 

1.0214167 

0.0018223 


2456284.54687 

-0.3960 

1.0129543 

0.0051269 

1.0212653 

-0.0083110 


2456284.54807 

-0.3933 

1.0188872 

0.0050400 

1.0211184 

-0.0022312 


2013 November 

2456605.75391 

-0.0035 

0.9584972 

0.0040311 

0.9926078 

-0.0341105 

-0.0272040 

2456605.75483 

-0.0014 

0.9682513 

0.0047093 

0.9919295 

-0.0236782 

-0.0273680 

2456605.75546 

0.0000 

0.9703735 

0.0040732 

0.9914676 

-0.0210941 

-0.0274150 

2456605.75639 

0.0021 

0.9600239 

0.0046811 

0.9907754 

-0.0307514 

-0.0273890 


Tololo Inter-American Observatory in Chile (CTIO), Sid¬ 
ing Spring Observatory in Australia (SSO), the South 
African Astronomical Observatory at Sutherland, South 
Africa (SAAO), and McDonald Observatory in Texas 
(MDO). The December 2012 data were obtained in sup¬ 
port of the Keck NIRSPEC observations and were in¬ 
tended to be simultaneous with the Keck Observations, 
but poor weather prevented data from being obtained 
on the exact night of observations. Instead, contempo¬ 
raneous data were obtained on 23 December 2012 UT. 
Exposure times were 90 sec, resulting in a cadence of 
approximately 103 seconds. All exposures were taken in 
the r' filter, similar to the filter used in the PTE-Orion 
survey data. 


After IBarnes et al.1 (|2013ll published the precession 
model, it was recognized that the orbital plane could 
process causing the transits to change, disappear and 
reappear, LCOGT was utilized again on 09 November 
2013 UT to observe PTFO 8-8695 data in order to con¬ 
firm the reappearance of the transit events and, in con¬ 
junction with the Spitzer photometry and the Keck ra¬ 
dial velocity data, confirm the precession model and the 
planetary nature of PTEO 8-8695b. The November 2013 
data were also obtained in the r' filter with an exposure 
time of 120 seconds; because the network of telescopes 
was able to observe the target with multiple telescopes 
at the same time, the effective cadence of the observa¬ 
tions ranged from w2 seconds to Ril30 seconds with the 
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BJD - 2456605.0 [days] 



Fig. 4. — Top: November 2013 LCOGT r' light curve, obtained approximately 3 years after the PTF data and 1.5 years after the Spitzer 
data, is plotted as a function of time. The blue circles indicate the data predicted to be within the transit window from the Spitzer 
ephemeris (Table [4||; the red points are outside the transit window and were used to fit a 8*^-order polynomial (black line) to model the 
stellar variability. Bottom: The LCOGT light curve, with the polynomial fit subtracted to remove the stellar variability, is plotted as a 
function of orbital phase. The solid black line represents the transit fit to the light curve and the data have been phased on this new 
ephemeris (see Table HJ. 


majority (« 50%) of the observations obtained at a 130 
second cadence. 

In order to create the differential photometry, an ini¬ 
tial normalization curve was created by taking the simple 
mean in magnitude space of all the raw reference light 
curves. For each reference star light curve, the mode of 
the residuals against the normalization curve was then 
subtracted, so that all the reference star light curves were 
normalized to the same flux level. Exposure by expo¬ 
sure, the mode of all the now-normalized reference star 
magnitudes was found, yielding the differential offset cor¬ 
rection needed for each exposure. These differential cor¬ 
rections were subtracted from the raw target light curve 
to produce the final differentially-corrected photometry. 
The same corrections applied to the original reference 
stars themselves (which should yield flat light curves) 
provided an internal consistency check. For a more de- 
tail ed overview of th e differ ential photometry technique, 
see Ivan Evken et al.l (|2011fl . The LCOGT light curves 
are tabulated in Tableland shown in Figures [3] and 01 

3. DISCUSSION 
3.1. Transit Ephemerides 

The Spitzer light curve was acquired 502 days past 
the time of the 2010 PTF transit midpoint determina¬ 
tion; in that time, the planet candidate PTFO 8-8695b 
would have orbited its host star 1119 times. As a re¬ 
sult, the timing of the transit in the Spitzer data was 


uncertain by 65 minutes. The transit was identified by 
visual inspection of the light curve at the time predicted 
by the PTF ephemeris. The transit event was found to 
be within the uncertainties of the predictions, but to lo¬ 
cate the transit more precisely, we excluded those data 
within ~2.5 hours of the predicted transit time (Fig. [1] 
blue circles) and fit an 8‘^’'-order polynomial to the data 
ontside the transit window (Fig. [T] red squares). Given 
the large number of data points in the Spitzer light curve 
outside the transit window (1174) and the fact that the 
polynomial is only intended to parametrically model the 
stellar variability, the order of the polynomial was set by 
requiring a non-significant change in the normalized vari¬ 
ance of the residuals. The polynomial fit was subtracted 
from the Spitzer light curve (data in the transit window 
included) to produce a whitened light curve suitable for 
a transit model fit (bottom of FigH]). 

The Spitzer light curve was fit using EXOFAST 
([Eastman et al.l [201311 which yielded a value and uncer¬ 
tainties via MCMC modeling for the Spitzer transit mid¬ 
point and the transit depth (see Table 0} . The Spitzer 
transit was found to be 17 minutes earlier than predicted, 
well within the 65 minute uncertainty from the 2010 PTF 
ephemeris uncertainty propagated to the time of the 2012 
Spitzer observations. Thus, more than 16 months after 
the PTF observations, the transit event occurred close 
to the predicted time. However, the transit, unlike the 
optical PTF transits, was found to be only ~ 0.6% deep. 
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TABLE 4 

Transit Ephemerides 



PTpa 

Spitzer 

LCOGT 

Approximate Date of Data Collection 

2009 Dec & 2010 Dec 

2012 April 

2013 November 

Transit Midpoint [Reduced BJD]^ 

5543.9402 

6045.7026 

6605.7555 

Midpoint Uncertainty 

±0.0008 

±0.0009 

±0.0004 

Transit Depth 

3-5% 

0.58 ± 0.01% 

2.74 ± 0.01% 

Days Past PTF Midpoint 


501.7622 

1061.8149 

Number of Orbits since PTF Midpoint'^ 


1119 

2368 

Uncertainty in Transit Prediction from PTF Midpoint 


±65 min. 

±136 min. 

Measured Offset from PTF Midpoint Prediction 


-17.2 min. 

—38.9 min. 

Days Past Spitzer Midpoint 



560.0527 

Number of Orbits since Spitzer Midpoint‘s 



1249 

Uncertainty in Transit Prediction from Spitzer Midpoint 



±71 min. 

Measured Offset from Spitzer Midpoint Prediction 



-21.8 min. 


^ From Ivan Evken et alj COia l ; IBarnes et al.l II2013I ') 

^ Reduced BJD = BJD - 2450000 
Period held constant at P = 0.448413 ± 0.000040 days. 


One of the primary purposes of the Spitzer observa¬ 
tions was to search for a secondary eclipse; however, 
we found no evidence for a secondary eclipse to within 
the limits of the data. At a semi-major axis distance of 
0.00838 AU from the star {Te/f = 3470K, = 1.4 i?©), 

the planet candidate should have an equilibrium temper¬ 
ature of Teq ~ 800 — 1800K, depending upon the albedo 
of the planet - similar to the expected effective tem- 

J ierature of a few m illion year old Jupiter-mass planet 
Baraffe et al.l l2003f l. The Spitzer 4.5pm light curve 
places a limit on the depth of the secondary eclipse of 
< 0.3% (3cr). This translates to an upper limit on the 
planetary candidate radius of Rp 1.5 Rjup, in reason¬ 
able agreement with the radius derived directly from the 
transit depths. 

The 2012 December LCOGT optical light curve was 
obtained 8 months after the Spitzer light curve, in sup¬ 
port of the Keck radial velocity data. In that time, 
the candidate planet PTFO 8-8695b would have orbited 
an additional 530 times and have a predicted transit 
time uncertainty of 10 minutes (0.02 in phase). As 
with the Spitzer light curve, data within ~2.5 hours of 
the predicted transit midpoint were excluded, and the 
out-of-transit data points were fit with a cubic spline; 
the flux brightening event that resembles a flare at 
BJD=2456284.625 prevented the use of a polynomial to 
parameterize the out-of-transit variability (Fig. [3]). After 
subtraction of the spline fit, the light curve shows no clear 
sign of a transit event to a limit of ^ 0.7% (Icr) within 
0.02 in phase of the predicted transit time. Structure is 
seen throughout the residual light curve at all phases and 
is likely the result of correlated noise in the data and/or 
true variability in the star. 

The 2013 November LCOGT optical light curve was 
obtained 560 days after the Spitzer light curve, as part of 
an effort to re-detect the transits - a necessary link in our 
efforts to confirm the orbital precession model. In that 
time, the the planetary candidate would have orbited an 
additional 1249 times and have a predicted transit time 
uncertainty of about 71 minutes. Even without the for¬ 
mal searching, the transit event is clearly seen at the time 
of predicted transit (see Fig. SI). To be consistent with 
our analysis of the Spitzer light curve and the previous 
LCOGT light curve, data within ^2.5 hours of the pre¬ 
dicted transit window were excluded (Fig. SI blue circles) 


and an 8*^-order polynomial was fit to the data outside 
the transit window (Fig. SI red squares). The polyno¬ 
mial fit was subtracted from the LCOGT light curve 
(including data within the transit window) to produce 
a whitened light curve suitable for a transit model fit 
(bottom of FigSl). 

The light curve was fit using EXOFAST which yielded 
a value for the LCOGT transit midpoint and the transit 
depth (see Table SJ- The transit was found 21 minutes 
earlier than predicted but again well within the 71 minute 
uncertainty from the Spitzer ephemeris uncertainty prop¬ 
agated to the time of the LCOGT observations; in fact, 
the measured time of transit was only 39 minutes ear¬ 
lier than what was predicted from the PTF ephemeris 
which was established 3 years and 2368 orbits earlier. 
Thus, the planetary candidate PTFO 8-8685b transited 
the star as predicted via simple extension of the transit 
timing ephemerides; this transit has a measured depth 
of ~ 2.7%. 

As a final note on the stability of the transit ephemeris, 
the p eriod was originally det ermined to within 3.5 sec¬ 
onds (jvan Evken et al.l 12012( 1 . If the measured transit 
midpoints for the PTF, Spitzer, and LCOGT light curves 
were perfect, the period uncertainty, propagated to the 
dates of the observations, would still produce an uncer¬ 
tainty greater than the measured offsets of the observed 
Spitzer and LCOGT transits. If the period was shorter 
by only 0.98 seconds, the transits would have been pre¬ 
dicted at the times observed, indicating the quality of 
the original measured ephemerides. 

3.2. Radial Velocity Limits 

When originally proposed, the 2012 December Keck 
NIRSPEC radial velocity observations were intended to 
measure the Rossiter-McLaughlin effect as the planet 
transited the star. The observations were timed to 
be centered on the transit midpoint, as predicted from 
the PTF ephemeris. The amplitude and shape of the 
Rossiter-McLaughlin effect depends heavily on the v sin i 
of the star, the size of the transiting planet, the orbital 
geometry of the system, and the viewing geometry from 
the Earth. We predicted that given the transit depth 
of ^ 4% and a stellar usini k, 80 km/s, the amplitude 
of the Rossiter-McLaughlin effect should be 2 — 4 km/s 
above the nominal radial velocity curve. 






















From Figured! it is clear that the Rossiter-McLaughlin 
effect was not detected, nor was the radial velocity signa¬ 
ture of the planetary candidate PTFO 8-8695b detected 
within the limits of the Keck NIRSPEC data. The radial 
velocity curve is consistent with a constant value with a 
reduced chi-squared = 0.85. The lack of detection 
of the Rossiter-McLaughlin effect is consistent with the 
LCOGT light curve from 2012 December which showed 
that there was not a transiting event at the time of the 
Keck observations. 

We have fitted the radial velocities with a Keplerian or¬ 
bit, where the orbital period was fixed to 0.448413 days 
from Ivan Evken et al.l (I2012D . and the orbital eccentric¬ 
ity wa s set to 0.0 consist e nt wi th both Ivan Evken et al.l 
(|2012f) and iBarnes et al.l (|2013f) . The time of the infe¬ 
rior conjunction (i.e., transit) was fixed by the Spitzer 
ephemeris listed in Table |4| The Keplerian orbital solu¬ 
tion to the radial velocity data in Table d) used RVLIN, 
a partially linearized, least-sq u ares f itting procedure de¬ 
scribed in iWright fc Howard (|2009ll : parameter uncer¬ 
tainties were estimated using the BOOTTRAN boot¬ 
strapping routines described in IWang et al.l (|2012f) . The 
results of the fitting are plotted over top the observations 
in Eig. di The Keplerian orbital fit does not improve sig¬ 
nificantly the fit over the constant weighted mean with 
only a reduced chi-squared of x^ = 0.80, and the semi¬ 
amplitude of the fit is statistically consistent with a flat 
line (K = 0.37 ± 0.33 km/s). 

3.3. Comparison to the Precession Models 

The mass of the host star PTFO 8-8695 was estimate d 
originally from isochrone fitting (|Briceho et al.l l2005f) . 
but the data c ould n ot di stinguish between models from 
iBaraffe et al.1 (j2010l l and iSiess et al.l (j200fll l: thus, leav- 
i ng an ambiguity in the stellar mass o f either ^ 0.44 
(jBaraffe et al l 1201011 or ^ 0.34 Mq (|Siess et al.l 120001 1. 
As a result. IBarnes et al.l 1)20131) developed a set of mod¬ 
els for both stellar mass estimates and fit the observed 
PTF transits from 2009 and 2010 separately and simul¬ 
taneously. The separate models described the expected 
rotational axis and orbital planet precession of the star 
and planet, and predicted the depth, shape, and times of 
the planetary transits for each of the stellar masses. Ad¬ 
ditionally, the precession modeling does take into account 
the wavelength dependence of the transits - in particu¬ 
lar, the effects of weakened limb and gravit y darkening at 
the w avelengths of the Spitzer observations (IBarnes et al.l 
I2013D . The data presented here a re compared to the two 
models from IBarnes et aP ()2013D . 

The higher stellar mass (0.44 Mq) precession model 
predicts a precession period of 581 days with a 184 day 
stretch with no transits visible at all as a result of the 
orbit processing out of the line of sight. Within the time 
that transits are visible, there is a double peak of deep 
transits separated by about 200 days. The lower mass 
(0.34 Mq) precession model predicts a similar transit 
structure as the higher mass model but on a more con¬ 
densed time scale (Figure [6]). The precession period is 
approximately 293 days and the time between the deep 
transits is about 80 days and the transits only disappear 
completely for approximately 75 days. 

Overall, the observations do not agree well with the 
0.44 Mq model. The Spitzer observations, observed to 
be only ^ 0.6% in depth are predicted by the model to 


be ^ 4%. Further, the observations in 2012 December 
are predicted to have occurred when the transit depth 
was at its deepest ^ 5%. Such a deep 5% deep tran¬ 
sit would have been detected in the LCOGT data which 
have a transit detection limit of ^ 0.7% (Fig. |3|). Addi¬ 
tionally, with a line-of-sight inclination oi i = 102°, the 
amplitude of the Rossiter-McLaughlin effect should have 
been > 2.5 km/s and would likely have been detected 
by the NIRSPEC radial velocity data at > 2a. Einally, 
the 2013 November optical transit detection is predicted 
at approximately the right depth, but the shape of the 
observed transit is much narrower in comparison to the 
predicted transit, indicating that the model incorrectly 
represents the system orientation at the time of the ob¬ 
servations. Thus, the 0.44 Mq precession model is incon¬ 
sistent with the observations for each of the three epochs 
of data presented. 

The 0.34 Mq precession model somewhat more closely 
matches the observations (Fig. Uj). The Spitzer light 
curve is represented well by the precession model transit 
prediction. Both the depth and the shape of the observed 
transit are in good agreement with the model predic¬ 
tion. Further, at the epoch of observations for the Keck 
NIRSPEC data and the first LCOGT light curve (2012 
December), the model predicts that the system should 
be non-transiting, and, indeed, the observations support 
the disappearance of the transiting events. The model 
predicts the return of the transits in late 2013, and the 
2013 November data confirm that the transits do reap¬ 
pear. However, the observations showed that the transits 
reappeared about 1 month earlier than predicted. With a 
precession orbital period of 293 days, the model is out-of- 
sync with the observations by ~ 10—15% of the predicted 
precession period. 

The cause of the inconsistency of the 0.34 Mq model 
with the 2013 November observations, while predicting 
the 2012 April and 2012 December observations correctly, 
is difficult to discern, but could be the result of a vari¬ 
ety of factors. The precession model depends heavily 
on knowledge of the stellar mass, the stellar moment of 
inertia, the stellar radius, and the planetary orbit eccen¬ 
tricity. The precession model assumes a circular orbit 
with a stellar radius of ~ 1.04 Rq] however, an eccentric 
orbit would allow for a larger stellar radius, which would 
affect the details of the precession model and the pre¬ 
cession period. The precession period is, to first order, 
proportional to the square of the stellar radius (i?*), and 
is also dependen t on higher order powers of the stellar 
radius (e.g., R^) ()Murrav fc Dermo^l2000t IBarnes et al.l 
IMl); a 5 — 7% change in the stellar radius could result 
in a 10 — 15% change in the precession period. 

The precession model also uses as inputs the orbital 
period and epochs of transits; even small errors on these, 
when propagated across three years and thousands of or¬ 
bits, may explain why the 0.34 Mq model predicts cor¬ 
rectly the 2012 data but not the 2013 data. Einally, the 
precession modeling fits not only the timing of the tran¬ 
sits but also the shape of the transits. If spots have sig¬ 
nificantly altered the shape of the transit, the precession 
model may not correctly predict the system parameters - 
particularly when the model is propagated through thou¬ 
sands of orbits. 

We cannot place better constraints on the stellar or 
planetary candidate parameters with the data presented 
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Fig. 5.— The precession model from [Harnes et all (1201311 for a stellar mass of M* = 0.44 Mq is compared to the observations. The 
bottom panel shows the predicted transit depths as a function of time starting with the 2009 PTF season and projected out through the 
end of 2014. The middle panel shows the transit model predictions at each of the times of observation. The PTF 2009 and PTF 2010 
transits match the model by default as these transits were used as constraints in the models. The transit model predictions at the time of 
the Spitzer, Keck/LCOGT, and LCOGT observations are shown. The top panel shows the orientation of the star and the planet during 
the time of observation. The shading on the star represents the brightness gradient caused by the gravity darkening; white is brighter and 
represents the poles of the star. The line-of-sight inclination (i) at the time of the observati on is also listed. Inc linations greater than 90° 
are related to the orientation of the system as a whole (for a description of the geometry see lBarnes et aOl2Q13l1 . 






Fig. 6 .— The precession model from IBarnes et al] 11201311 for a stellar mass of M* = 0.34 Mq is compared to the observations. The figure 
content is the same as for Figure [5] 
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here. The parameter space of the precession models is 
extremely large and is difficult to narrow and beyond the 
intended scope of this paper. A recent paper, studying 
the same system with our PTF data but without requir¬ 
ing the spin-orbit to be synchronously locked, find that 
a precession model can still reproduce the observations, 
with a precessi on period of 199 days and a planet mass of 
^ 4 — 5Mjup (|Kamiaka et al .11^01511 . We are continuing 
to observe the planetary candidate with long-term pho¬ 
tometric observations and radial velocity monitoring in 
an effort to assess more fully the existence and planetary 
nature of PTFO 8-8695b and to more carefully determine 
the validity and parameters of the precession model. 

While the observations presented here cannot uniquely 
identify which precession model might best describe the 
data, we can utilize the NIRSPEC radial velocity up¬ 
per limits on the semi-amplitude to set limits on the 
potential mass of the planetary candidate. If the stel¬ 
lar mass is 0.44 Mq, then the mass limit on a possi¬ 
ble planetary candidate is Mp sini < 0.8 Mjup] if the 
stellar mass is closer to 0.34 Mq, then the mass limit 
is Mp sin i < 0.7 Mjup. As indicated by the data, the 
planet candidate does not transit during the time of the 
radial velocity observations, which sets an upper limit on 
the line-of-sight orbital inclination of i < 50°. At that 
inclination, the planet candidate mass would have a mass 
of Mp < 0.9 — 1.1 Mjup. If the orbital inclination is as 
low as i ^ 20° — 30° as predicted by the 0.34 Mq pre¬ 
cession model, then the planetary candidate could have 
a mass of Mp < 2.0 — 2.5 Mjup. 

3.4. Stellar Activity 

PTFO 8-8685 is a low mass young star and as a re- 
sult is variable and is likely s potted as d iscussed in 
Ivan Evken et al.l (|2012D fsee also iKoeiil 120151) . The op¬ 
tical radial velocity curve modulations were attributed 
to the presence of spots, although not necessarily to the 
presence of a single spot. Large and long-lived, high 
latitude spots ar e commonly presen t in weak-l i ned T 
Tauri stars (e.g., iStelzer et al.l 1200^ iRice et rdl 1201111 . 
and spots have been seen to mimic eclipsi ng companions 
in ot her systems (e.g., R,X,11608.6-3922: Lloergens et al.l 
1200111 . In that work, light curve data obtained in multiple 
filters and spectroscopic monitoring enabled the authors 
to determine that the deep eclipses (0.5 mag), with long 
eclipse durations (^ 0.5 in phase), were caused by stellar 
spots which disappeared in observations obtained 4 years 
later. They also found that the photometric variations 
were also consistent with the 2.4 km/s radial velocity 
v ariations but not compa tible with an eclipsing binary. 

Ivan Evken et al.l (|2012ll could not find a spot distribu¬ 
tion solution to ht the observed light curves for PTFO 8- 
8695, but we revisit the possibility of spots as a cause 
of the observed transit events here. Unlike RXJ1608, 
the observed transit features for PTFO 8-8695 are more 
shallow (3% —5%) and shorter in duration (0.1 in phase), 
and, as seen in 2009, sometimes flat-bottom during the 
event. In principle, a spotted surface with a single hot 
spot near the pole briefly eclipsed by the limb of the 
star may be able produce the short transit duration ob¬ 
served. However, such a hot spot would be expected 
to show a brightening 180° out of phase with the flux 
dip but such a brightening is not seen in the data. In 
fact, many of the light curves from the 2009/2010 dis¬ 


covery data (jvan Evken et al.l 1 20 121) display a brightness 
decrease near these phases in the light curves. 

A cold spot below the equator could briefly come into 
view. In general, the spot size (i.e., the surface filling- 
factor) needed to produce the observed transit depth can 
be estimated if we assume a photospheric temperature to 
spot temperature ratio {Tp/Tg). T he star has an effec- 
tive temperature of Tg// ~ 3470K (|Briceno et al.l[2005[l 
and M-dwarf spots can be 500 — lOOOK coo ler than the 
photosphere (| Barnes. Jeffers, fc .lone^l^Ollh . The opti¬ 
cal (0.65 /rm) detections of the transit event have depths 
that range from 3% — 5% corresponding to spot filling 
factors of ^ 3.5% — 8%. After taking into account the 
reduced contrast with the photosphere at longer wave¬ 
lengths, the 0.58% infrared (4.5/im) transit depth cor¬ 
responds to a filling-factor coverage of 1.5% — 2.5% 
- a factor of 2 — 4 times smaller. To explain the short 
transition duration of 0.1 in phase (^ 36° in longitude), 
such a spot would need to have a viewable longitudinal 
extent of ~ 36/2 ~ 18°, and to yield a spot that covers 
2% — 8% of the surface, the viewable latitudinal extent 
of the spot would need to be 20° — 40°, but a spot of 
such latitudinal extent would likely be visible for more 
than just 1/10 of the rotation period. 

All of this is not impossible for an active low-mass star, 
but the transit events would need to change shape and 
depth (and be flat-bottomed as is the case for the 2009 
data). At the same time, the spots would need to disap¬ 
pear and reappear at nearly the same stellar longitude 
over the course of > 3 years. As indicated in section ITTI 
the transit times agree with the predicted ephemerides 
to within < 20 minutes or about < 3% of the period. 
The spot would need to appear, evolve, disappear, and 
re-appear all within < 10° of the same stellar longitudi¬ 
nal position as the event timing ephemerides have been 
consistent for over three years, and do this all within the 
context of a generally spotted star producing the contin¬ 
uous va riability of the sta r. _ 

As in Ivan Evken et al.l (I2012D . we find it difficult to 
model the transit events in a self-consistent manner with 
the data acquired over 3 years. We are, however, contin¬ 
uing to pursue multi-year and multi-color light curves in 
an effort to understand better the nature of the stellar 
variability. 

In section 13.11 we refer to the brightening event in 
the 2012 December optical light curve that occurs at 
phase 0.8. This brightening does not look like a typ- 
ical flare with a shar p rise and exponential decline 
l|Walkowicz et al.l 1201111 . One possibility is that this 
brightening is the result of an accretion event on the 
star. PTFO 8-8695 exhibits relatively strong H-a emis¬ 
sion and it has been noted that the planetary candidate 
is near or at the Roche- limit - particularly, if the pu¬ 
tative orbit is eccentric (|van Evken et al.ll2012fl . As a 
result, the planetary candidate may be evaporating and 
the “flare” may actually be the result of infall onto the 
stellar surface. More detailed work on the H-a emis- 
sion and its variability is the subject of another paper 
(j.Iohns-Krull et aDl20i5[l . 

4. SUMMARY 

PTFO 8-8695b was discovered in the PTF-Orion sur¬ 
vey for transiting exoplanets bv Ivan Evken et ^ (|2012ll : 
PTFO 8-8695b was found to be in an 0.45 day orbit 
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and to have a mass of ;< 4 — 5M j,ir,. That discovery 
was followed-up by a prediction by iBarnes et al.l (|2013[i 
that the stellar rotation axis and the orbital plane of the 
planet should process and that the transits should change 
shape and depth and disappear and reappear with a pe¬ 
riod of 300 — 600 days. The two papers put together a 
coherent picture of a Jupiter-sized planetary candidate 
orbiting a ^3MYr old weak-lined T Tauri star; however, 
the precession model needed confirmation and PTFO 8- 
8695b remained a planetary candidate with a mass less 
than < ZMjup- 

We have obtained follow-up observations of PTFO 8- 
8695 that includes Spitzer 4.5/rm light curve observa¬ 
tions, Keck NIRSPEC radial velocity observations, and 
LCOGT optical light curve observations. The data con¬ 
firm that the transit events are persistent over many 
years and the times of the transits are predictable 
from the transit timing measurements, consistent with 
a steady period of 0.448 days. The transit events do ap¬ 
pear to grow more shallow, disappear, and reappear as 
predicted by the precession modeling. 

However, the precession model and observations are 
not in perfect agreement and more observations are 
needed to place better constraints on the model and to 
confirm the planetary nature of the candidate PTFO 8- 
9695b. Towards this end, we are pursuing additional 
long-term transit observations with LCOGT and radial 
velocity monitoring that will help limit the models, the 
orbit, the stellar parameters, and help confirm or refute 
the planetary nature as the source of the observed tran¬ 
sit events which have remained for over three years of 
observations. 
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